Autosomal dominant in cis D4Z4 repeat array duplication alleles in facioscapulohumeral dystrophy.

Facioscapulohumeral dystrophy (FSHD) has a unique genetic aetiology resulting in partial chromatin relaxation of the D4Z4 macrosatellite repeat array on 4qter. This D4Z4 chromatin relaxation facilitates inappropriate expression of the transcription factor DUX4 in skeletal muscle. DUX4 is encoded by a retrogene that is embedded within the distal region of the D4Z4 repeat array. In the European population, the D4Z4 repeat array is usually organized in a single array that ranges between 8 and 100 units. D4Z4 chromatin relaxation and DUX4 derepression in FSHD is most often caused by repeat array contraction to 1-10 units (FSHD1) or by a digenic mechanism requiring pathogenic variants in a D4Z4 chromatin repressor like SMCHD1, combined with a repeat array between 8 and 20 units (FSHD2). With a prevalence of 1.5% in the European population, in cis duplications of the D4Z4 repeat array, where two adjacent D4Z4 arrays are interrupted by a spacer sequence, are relatively common but their relationship to FSHD is not well understood. In cis duplication alleles were shown to be pathogenic in FSHD2 patients; however, there is inconsistent evidence for the necessity of an SMCHD1 mutation for disease development. To explore the pathogenic nature of these alleles we compared in cis duplication alleles in FSHD patients with or without pathogenic SMCHD1 variant. For both groups we showed duplication-allele-specific DUX4 expression. We studied these alleles in detail using pulsed-field gel electrophoresis-based Southern blotting and molecular combing, emphasizing the challenges in the characterization of these rearrangements. Nanopore sequencing was instrumental to study the composition and methylation of the duplicated D4Z4 repeat arrays and to identify the breakpoints and the spacer sequence between the arrays. By comparing the composition of the D4Z4 repeat array of in cis duplication alleles in both groups, we found that specific combinations of proximal and distal repeat array sizes determine their pathogenicity. Supported by our algorithm to predict pathogenicity, diagnostic laboratories should now be furnished to accurately interpret these in cis D4Z4 repeat array duplications, alleles that can easily be missed in routine settings.

Deciphering D4Z4 CpG methylation gradients in fascioscapulohumeral muscular dystrophy using nanopore sequencing.

Fascioscapulohumeral muscular dystrophy (FSHD) is caused by a unique genetic mechanism that relies on contraction and hypomethylation of the D4Z4 macrosatellite array on the Chromosome 4q telomere allowing ectopic expression of the DUX4 gene in skeletal muscle. Genetic analysis is difficult because of the large size and repetitive nature of the array, a nearly identical array on the 10q telomere, and the presence of divergent D4Z4 arrays scattered throughout the genome. Here, we combine nanopore long-read sequencing with Cas9-targeted enrichment of 4q and 10q D4Z4 arrays for comprehensive genetic analysis including determination of the length of the 4q and 10q D4Z4 arrays with base-pair resolution. In the same assay, we differentiate 4q from 10q telomeric sequences, determine A/B haplotype, identify paralogous D4Z4 sequences elsewhere in the genome, and estimate methylation for all CpGs in the array. Asymmetric, length-dependent methylation gradients were observed in the 4q and 10q D4Z4 arrays that reach a hypermethylation point at approximately 10 D4Z4 repeat units, consistent with the known threshold of pathogenic D4Z4 contractions. High resolution analysis of individual D4Z4 repeat methylation revealed areas of low methylation near the CTCF/insulator region and areas of high methylation immediately preceding the DUX4 transcriptional start site. Within the DUX4 exons, we observed a waxing/waning methylation pattern with a 180-nucleotide periodicity, consistent with phased nucleosomes. Targeted nanopore sequencing complements recently developed molecular combing and optical mapping approaches to genetic analysis for FSHD by adding precision of the length measurement, base-pair resolution sequencing, and quantitative methylation analysis.

Promising AAV.U7snRNAs vectors targeting DMPK improve DM1 hallmarks in patient-derived cell lines.

Myotonic dystrophy type 1 (DM1) is the most common form of muscular dystrophy in adults and affects mainly the skeletal muscle, heart, and brain. DM1 is caused by a CTG repeat expansion in the 3'UTR region of the DMPK gene that sequesters muscleblind-like proteins, blocking their splicing activity and forming nuclear RNA foci. Consequently, many genes have their splicing reversed to a fetal pattern. There is no treatment for DM1, but several approaches have been explored, including antisense oligonucleotides (ASOs) aiming to knock down DMPK expression or bind to the CTGs expansion. ASOs were shown to reduce RNA foci and restore the splicing pattern. However, ASOs have several limitations and although being safe treated DM1 patients did not demonstrate improvement in a human clinical trial. AAV-based gene therapies have the potential to overcome such limitations, providing longer and more stable expression of antisense sequences. In the present study, we designed different antisense sequences targeting exons 5 or 8 of DMPK and the CTG repeat tract aiming to knock down DMPK expression or promote steric hindrance, respectively. The antisense sequences were inserted in U7snRNAs, which were then vectorized in AAV8 particles. Patient-derived myoblasts treated with AAV8. U7snRNAs showed a significant reduction in the number of RNA foci and re-localization of muscle-blind protein. RNA-seq analysis revealed a global splicing correction in different patient-cell lines, without alteration in DMPK expression.

A genome-wide association analysis of loss of ambulation in dystrophinopathy patients suggests multiple candidate modifiers of disease severity.

The major determinant of disease severity in Duchenne muscular dystrophy (DMD) or milder Becker muscular dystrophy (BMD) is whether the dystrophin gene (DMD) mutation truncates the mRNA reading frame or allows expression of a partially functional protein. However, even in the complete absence of dystrophin, variability in disease severity is observed, and candidate gene studies have implicated several genes as modifiers. Here we present the largest genome-wide search to date for loci influencing severity in N = 419 DMD patients. Availability of subjects for such studies is quite limited, leading to modest sample sizes, which present a challenge for GWAS design. We have therefore taken special steps to minimize heterogeneity within our dataset at the DMD locus itself, taking a novel approach to mutation classification to effectively exclude the possibility of residual dystrophin expression, and utilized statistical methods that are well adapted to smaller sample sizes, including the use of a novel linear regression-like residual for time to ambulatory loss and the application of evidential statistics for the GWAS approach. Finally, we applied an unbiased in silico pipeline, utilizing functional genomic datasets to explore the potential impact of the best supported SNPs. In all, we obtained eight SNPs (out of 1,385,356 total) with posterior probability of trait-marker association (PPLD) ≥ 0.4, representing six distinct loci. Our analysis prioritized likely non-coding SNP regulatory effects on six genes (ETAA1, PARD6G, GALNTL6, MAN1A1, ADAMTS19, and NCALD), each with plausibility as a DMD modifier. These results support both recurrent and potentially new pathways for intervention in the dystrophinopathies.

Deciphering D4Z4 CpG methylation gradients in fascioscapulohumeral muscular dystrophy using nanopore sequencing.

Fascioscapulohumeral muscular dystrophy (FSHD) is caused by a unique genetic mechanism that relies on contraction and hypomethylation of the D4Z4 macrosatellite array on the chromosome 4q telomere allowing ectopic expression of the DUX4 gene in skeletal muscle. Genetic analysis is difficult due to the large size and repetitive nature of the array, a nearly identical array on the 10q telomere, and the presence of divergent D4Z4 arrays scattered throughout the genome. Here, we combine nanopore long-read sequencing with Cas9-targeted enrichment of 4q and 10q D4Z4 arrays for comprehensive genetic analysis including determination of the length of the 4q and 10q D4Z4 arrays with base-pair resolution. In the same assay, we differentiate 4q from 10q telomeric sequences, determine A/B haplotype, identify paralogous D4Z4 sequences elsewhere in the genome, and estimate methylation for all CpGs in the array. Asymmetric, length-dependent methylation gradients were observed in the 4q and 10q D4Z4 arrays that reach a hypermethylation point at approximately 10 D4Z4 repeat units, consistent with the known threshold of pathogenic D4Z4 contractions. High resolution analysis of individual D4Z4 repeat methylation revealed areas of low methylation near the CTCF/insulator region and areas of high methylation immediately preceding the DUX4 transcriptional start site. Within the DUX4 exons, we observed a waxing/waning methylation pattern with a 180-nucleotide periodicity, consistent with phased nucleosomes. Targeted nanopore sequencing complements recently developed molecular combing and optical mapping approaches to genetic analysis for FSHD by adding precision of the length measurement, base-pair resolution sequencing, and quantitative methylation analysis.

Systemic delivery of an AAV9 exon-skipping vector significantly improves or prevents features of Duchenne muscular dystrophy in the Dup2 mouse.

Duchenne muscular dystrophy (DMD) is typically caused by mutations that disrupt the DMD reading frame, but nonsense mutations in the 5' part of the gene induce utilization of an internal ribosomal entry site (IRES) in exon 5, driving expression of a highly functional N-truncated dystrophin. We have developed an AAV9 vector expressing U7 small nuclear RNAs targeting DMD exon 2 and have tested it in a mouse containing a duplication of exon 2, in which skipping of both exon 2 copies induces IRES-driven expression, and skipping of one copy leads to wild-type dystrophin expression. One-time intravascular injection either at postnatal days 0-1 or at 2 months results in efficient exon skipping and dystrophin expression, and significant protection from functional and pathologic deficits. Immunofluorescence quantification showed 33%-53% average dystrophin intensity and 55%-79% average dystrophin-positive fibers in mice treated in adulthood, with partial amelioration of DMD pathology and correction of DMD-associated alterations in gene expression. In mice treated neonatally, dystrophin immunofluorescence reached 49%-85% of normal intensity and 76%-99% dystrophin-positive fibers, with near-complete correction of dystrophic pathology, and these beneficial effects persisted for at least 6 months. Our results demonstrate the robustness, durability, and safety of exon 2 skipping using scAAV9.U7snRNA.ACCA, supporting its clinical use.

Intron mutations and early transcription termination in Duchenne and Becker muscular dystrophy.

DMD pathogenic variants for Duchenne and Becker muscular dystrophy are detectable with high sensitivity by standard clinical exome analyses of genomic DNA. However, up to 7% of DMD mutations are deep intronic and analysis of muscle-derived RNA is an important diagnostic step for patients who have negative genomic testing but abnormal dystrophin expression in muscle. In this study, muscle biopsies were evaluated from 19 patients with clinical features of a dystrophinopathy, but negative clinical DMD mutation analysis. Reverse transcription-polymerase chain reaction or high-throughput RNA sequencing methods identified 19 mutations with one of three pathogenic pseudoexon types: deep intronic point mutations, deletions or insertions, and translocations. In association with point mutations creating intronic splice acceptor sites, we observed the first examples of DMD pseudo 3'-terminal exon mutations causing high efficiency transcription termination within introns. This connection between splicing and premature transcription termination is reminiscent of U1 snRNP-mediating telescripting in sustaining RNA polymerase II elongation across large genes, such as DMD. We propose a novel classification of three distinct types of mutations identifiable by muscle RNA analysis, each of which differ in potential treatment approaches. Recognition and appropriate characterization may lead to therapies directed toward full-length dystrophin expression for some patients.

Phenotypic Spectrum of Dystrophinopathy Due to Duchenne Muscular Dystrophy Exon 2 Duplications.

BACKGROUND AND OBJECTIVES: To describe the phenotypic spectrum of dystrophinopathy in a large cohort of individuals with DMD exon 2 duplications (Dup2), who may be particularly amenable to therapies directed at restoring expression of either full-length dystrophin or nearly full-length dystrophin through utilization of the DMD exon 5 internal ribosome entry site (IRES). METHODS: In this retrospective observational study, we analyzed data from large genotype-phenotype databases (the United Dystrophinopathy Project [UDP] and the Italian DMD network) and classified participants into Duchenne muscular dystrophy (DMD), intermediate muscular dystrophy (IMD), or Becker muscular dystrophy (BMD) phenotypes. Log-rank tests for time-to-event variables were used to compare age at loss of ambulation (LOA) in participants with Dup2 vs controls without Dup2 in the UDP database and for comparisons between steroid-treated vs steroid-naive participants with Dup2. RESULTS: Among 66 participants with Dup2 (UDP = 40, Italy = 26), 61% were classified as DMD, 9% as IMD, and 30% as BMD. Median age at last observation was 15.4 years (interquartile range 8.79-26.0) and 75% had been on corticosteroids for at least 6 months. Age at LOA differed significantly between participants with Dup2 DMD and historical controls without Dup2 DMD (p < 0.001). Valid spirometry was limited but suggested a delay in the typical age-related decline in forced vital capacity and 24 of 55 participants with adequate cardiac data had cardiomyopathy. DISCUSSION: Some patients with Dup2 display a milder disease course than controls without Dup2 DMD, and prolonged ambulation with corticosteroids suggests the potential of IRES activation as a molecular mechanism. As Dup2-targeted therapies reach clinical applications, this information is critical to aid in the interpretation of the efficacy of new treatments.

Absence of Significant Off-Target Splicing Variation with a U7snRNA Vector Targeting DMD Exon 2 Duplications.

Exon skipping therapies for Duchenne muscular dystrophy that restore an open reading frame can be induced by the use of noncoding U7 small nuclear RNA (U7snRNA) modified by an antisense exon-targeting sequence delivered by an adeno-associated virus (AAV) vector. We have developed an AAV vector (AAV9.U7-ACCA) containing four U7snRNAs targeting the splice donor and acceptor sites of dystrophin exon 2, resulting in highly efficient exclusion of DMD exon 2. We assessed the specificity of splice variation induced by AAV9.U7-ACCA delivery in the Dmd exon 2 duplication (Dup2) mouse model through an unbiased RNA-seq approach. Treatment-related effects on pre-mRNA splicing were quantified using local splicing variation (LSV) analysis. Filtering the transcriptome for differences in treatment-related splicing resulted in only 16 candidate off-target LSVs. Only a single candidate off-target LSV was found in both skeletal and cardiac muscle tissue and occurred at a known variable cassette exon. In contrast, four LSVs represented significant on-target correction of Dmd exon 2 splicing and transcriptome analysis showed correction of known dystrophin-deficient gene dysregulation. We conclude that the absence of off-target splicing induced by treatment with the U7-ACCA vector supports the continued clinical development of this approach.

Direct Reprogramming of Human Fibroblasts into Myoblasts to Investigate Therapies for Neuromuscular Disorders.

Investigations into both the pathophysiology and therapeutic targets in muscular dystrophies have been hampered by the limited proliferative capacity of human myoblasts. Several mouse models have been created but they either do not truly represent the human physiopathology of the disease or are not representative of the broad spectrum of mutations found in humans. The immortalization of human primary myoblasts is an alternative to this limitation; however, it is still dependent on muscle biopsies, which are invasive and not easily available. In contrast, skin biopsies are easier to obtain and less invasive to patients. Fibroblasts derived from skin biopsies can be immortalized and transdifferentiated into myoblasts, providing a source of cells with excellent myogenic potential. Here, we describe a fast and direct reprogramming method of fibroblast into a myogenic lineage. Fibroblasts are transduced with two lentiviruses: hTERT to immortalize the primary culture and a tet-inducible MYOD, which upon the addition of doxycycline, induces the conversion of fibroblasts into myoblasts and then mature myotubes, which express late differentiation markers. This quick transdifferentiation protocol represents a powerful tool to investigate pathological mechanisms and to investigate innovative gene-based or pharmacological biotherapies for neuromuscular disorders.

Incipient genome erosion and metabolic streamlining for antibiotic production in a defensive symbiont.

Genome erosion is a frequently observed result of relaxed selection in insect nutritional symbionts, but it has rarely been studied in defensive mutualisms. Solitary beewolf wasps harbor an actinobacterial symbiont of the genus Streptomyces that provides protection to the developing offspring against pathogenic microorganisms. Here, we characterized the genomic architecture and functional gene content of this culturable symbiont using genomics, transcriptomics, and proteomics in combination with in vitro assays. Despite retaining a large linear chromosome (7.3 Mb), the wasp symbiont accumulated frameshift mutations in more than a third of its protein-coding genes, indicative of incipient genome erosion. Although many of the frameshifted genes were still expressed, the encoded proteins were not detected, indicating post-transcriptional regulation. Most pseudogenization events affected accessory genes, regulators, and transporters, but "Streptomyces philanthi" also experienced mutations in central metabolic pathways, resulting in auxotrophies for biotin, proline, and arginine that were confirmed experimentally in axenic culture. In contrast to the strong A+T bias in the genomes of most obligate symbionts, we observed a significant G+C enrichment in regions likely experiencing reduced selection. Differential expression analyses revealed that-compared to in vitro symbiont cultures-"S. philanthi" in beewolf antennae showed overexpression of genes for antibiotic biosynthesis, the uptake of host-provided nutrients and the metabolism of building blocks required for antibiotic production. Our results show unusual traits in the early stage of genome erosion in a defensive symbiont and suggest tight integration of host-symbiont metabolic pathways that effectively grants the host control over the antimicrobial activity of its bacterial partner.

High throughput screening for expanded CTG repeats in myotonic dystrophy type 1 using melt curve analysis.

BACKGROUND: Myotonic dystrophy type 1 (DM1) is caused by CTG repeat expansions in the DMPK gene and is the most common form of muscular dystrophy. Patients can have long delays from onset to diagnosis, since clinical signs and symptoms are often nonspecific and overlapping with other disorders. Clinical genetic testing by Southern blot or triplet-primed PCR (TP-PCR) is technically challenging and cost prohibitive for population surveys. METHODS: Here, we present a high throughput, low-cost screening tool for CTG repeat expansions using TP-PCR followed by high resolution melt curve analysis with saturating concentrations of SYBR GreenER dye. RESULTS: We determined that multimodal melt profiles from the TP-PCR assay are a proxy for amplicon length stoichiometry. In a screen of 10,097 newborn blood spots, melt profile analysis accurately reflected the tri-modal distribution of common alleles from 5 to 35 CTG repeats, and identified the premutation and full expansion alleles. CONCLUSION: We demonstrate that robust detection of expanded CTG repeats in a single tube can be achieved from samples derived from specimens with minimal template DNA such as dried blood spots (DBS). This technique is readily adaptable to large-scale testing programs such as population studies and newborn screening programs.

Core transcriptional networks in Williams syndrome: IGF1-PI3K-AKT-mTOR, MAPK and actin signaling at the synapse echo autism.

Gene networks for disorders of social behavior provide the mechanisms critical for identifying therapeutic targets and biomarkers. Large behavioral phenotypic effects of small human deletions make the positive sociality of Williams syndrome (WS) ideal for determining transcriptional networks for social dysfunction currently based on DNA variations for disorders such as autistic spectrum disorder (ASD) and schizophrenia (SCHZ). Consensus on WS networks has been elusive due to the need for larger cohort size, sensitive genome-wide detection and analytic tools. We report a core set of WS network perturbations in a cohort of 58 individuals (34 with typical, 6 atypical deletions and 18 controls). Genome-wide exon-level expression arrays robustly detected changes in differentially expressed gene (DEG) transcripts from WS deleted genes that ranked in the top 11 of 12 122 transcripts, validated by quantitative reverse transcription PCR, RNASeq and western blots. WS DEG's were strictly dosed in the full but not the atypical deletions that revealed a breakpoint position effect on non-deleted CLIP2, a caveat for current phenotypic mapping based on copy number variants. Network analyses tested the top WS DEG's role in the dendritic spine, employing GeneMANIA to harmonize WS DEGs with comparable query gene-sets. The results indicate perturbed actin cytoskeletal signaling analogous to the excitatory dendritic spines. Independent protein-protein interaction analyses of top WS DEGs generated a 100-node graph annotated topologically revealing three interacting pathways, MAPK, IGF1-PI3K-AKT-mTOR/insulin and actin signaling at the synapse. The results indicate striking similarity of WS transcriptional networks to genome-wide association study-based ASD and SCHZ risk suggesting common network dysfunction for these disorders of divergent sociality.

Population-Based Prevalence of Myotonic Dystrophy Type 1 Using Genetic Analysis of Statewide Blood Screening Program.

OBJECTIVE: To determine whether the genetic prevalence of the CTG expansion in the DMPK gene associated with myotonic dystrophy type 1 (DM1) in an unbiased cohort is higher than previously reported population estimates, ranging from 5 to 20 per 100,000 individuals. METHODS: This study used a cross-sectional cohort of deidentified dried blood spots from the newborn screening program in the state of New York, taken from consecutive births from 2013 to 2014. Blood spots were screened for the CTG repeat expansion in the DMPK gene using triplet-repeat primed PCR and melt curve analysis. Melt curve morphology was assessed by 4 blinded reviewers to identify samples with possible CTG expansion. Expansion of the CTG repeat was validated by PCR fragment sizing using capillary electrophoresis for samples classified as positive or premutation to confirm the result. Prevalence was calculated as the number of samples with CTG repeat size ≥50 repeats compared to the overall cohort. RESULTS: Of 50,382 consecutive births, there were 24 with a CTG repeat expansion ≥50, consistent with a diagnosis of DM1. This represents a significantly higher DM1 prevalence of 4.76 per 10,000 births (95% confidence interval 2.86-6.67) or 1 in every 2,100 births. There were an additional 96 samples (19.1 per 10,000 or 1 in 525 births) with a CTG expansion in the DMPK gene in the premutation range (CTG)35-49. CONCLUSION: The prevalence of individuals with CTG repeat expansions in DMPK is up to 5 times higher than previous reported estimates. This suggests that DM1, with multisystemic manifestations, is likely underdiagnosed in practice.

X-linked muscular dystrophy in a Labrador Retriever strain: phenotypic and molecular characterisation.

BACKGROUND: Canine models of Duchenne muscular dystrophy (DMD) are a valuable tool to evaluate potential therapies because they faithfully reproduce the human disease. Several cases of dystrophinopathies have been described in canines, but the Golden Retriever muscular dystrophy (GRMD) model remains the most used in preclinical studies. Here, we report a new spontaneous dystrophinopathy in a Labrador Retriever strain, named Labrador Retriever muscular dystrophy (LRMD). METHODS: A colony of LRMD dogs was established from spontaneous cases. Fourteen LRMD dogs were followed-up and compared to the GRMD standard using several functional tests. The disease causing mutation was studied by several molecular techniques and identified using RNA-sequencing. RESULTS: The main clinical features of the GRMD disease were found in LRMD dogs; the functional tests provided data roughly overlapping with those measured in GRMD dogs, with similar inter-individual heterogeneity. The LRMD causal mutation was shown to be a 2.2-Mb inversion disrupting the DMD gene within intron 20 and involving the TMEM47 gene. In skeletal muscle, the Dp71 isoform was ectopically expressed, probably as a consequence of the mutation. We found no evidence of polymorphism in either of the two described modifier genes LTBP4 and Jagged1. No differences were found in Pitpna mRNA expression levels that would explain the inter-individual variability. CONCLUSIONS: This study provides a full comparative description of a new spontaneous canine model of dystrophinopathy, found to be phenotypically equivalent to the GRMD model. We report a novel large DNA mutation within the DMD gene and provide evidence that LRMD is a relevant model to pinpoint additional DMD modifier genes.

Clinical Phenotypes of DMD Exon 51 Skip Equivalent Deletions: A Systematic Review.

BACKGROUND: Eteplirsen, the first FDA-approved RNA-modifying therapy for DMD, is applicable to ∼13% of patients with DMD. Because multiple exonic deletions are amenable to exon 51 skipping, the isoforms resulting from the various exon 51-skipped transcripts may vary in stability, function, and phenotype. OBJECTIVE/METHODS: We conducted a detailed review of dystrophinopathy published literature and unpublished databases to compile phenotypic features of patients with exon 51 "skip-equivalent" deletions. RESULTS: Theoretically, 48 different in-frame transcripts may result from exon 51 skipping. We found sufficient clinical information on 135 patients carrying mutations that would result in production of 11 (23%) of these transcripts, suggesting the remainder have not been identified in vivo. The majority had mild phenotypes: BMD (n = 81) or isolated dilated cardiomyopathy (n = 3). Particularly interesting are the asymptomatic (n = 10) or isolated hyperCKemia (n = 20) patients with deletions of exons 45- 51, 48- 51, 49- 51 and 50- 51. Finally, 16 (12%) had more severe phenotypes described as intermediate (n = 2) or DMD (n = 14), and 6 reports had no definitive phenotype. CONCLUSIONS: This review shows that the majority of exon 51 "skip-equivalent" deletions result in milder (BMD) phenotypes and supports that exon 51 skipping therapy could provide clinical benefit, although we acknowledge that other factors, such as age at treatment initiation or ongoing standard of care, may influence the degree of benefit.

Clinical, genetic, and pathologic characterization of FKRP Mexican founder mutation c.1387A>G.

OBJECTIVE: To characterize the clinical phenotype, genetic origin, and muscle pathology of patients with the FKRP c.1387A>G mutation. METHODS: Standardized clinical data were collected for all patients known to the authors with c.1387A>G mutations in FKRP. Muscle biopsies were reviewed and used for histopathology, immunostaining, Western blotting, and DNA extraction. Genetic analysis was performed on extracted DNA. RESULTS: We report the clinical phenotypes of 6 patients homozygous for the c.1387A>G mutation in FKRP. Onset of symptoms was <2 years, and 5 of the 6 patients never learned to walk. Brain MRIs were normal. Cognition was normal to mildly impaired. Microarray analysis of 5 homozygous FKRP c.1387A>G patients revealed a 500-kb region of shared homozygosity at 19q13.32, including FKRP. All 4 muscle biopsies available for review showed end-stage dystrophic pathology, near absence of glycosylated α-dystroglycan (α-DG) by immunofluorescence, and reduced molecular weight of α-DG compared with controls and patients with homozygous FKRP c.826C>A limb-girdle muscular dystrophy. CONCLUSIONS: The clinical features and muscle pathology in these newly reported patients homozygous for FKRP c.1387A>G confirm that this mutation causes congenital muscular dystrophy. The clinical severity might be explained by the greater reduction in α-DG glycosylation compared with that seen with the c.826C>A mutation. The shared region of homozygosity at 19q13.32 indicates that FKRP c.1387A>G is a founder mutation with an estimated age of 60 generations (∼1,200-1,500 years).

A genome-wide association study for mastitis resistance in phenotypically well-characterized Holstein dairy cattle using a selective genotyping approach.

A decrease in the incidence of bovine mastitis, the costliest disease in the dairy industry, can be facilitated through genetic marker-assisted selective breeding programs. Identification of genomic variants associated with mastitis resistance is an ongoing endeavor for which genome-wide association studies (GWAS) using high-density arrays provide a valuable tool. We identified single nucleotide polymorphisms (SNPs) in Holstein dairy cattle associated with mastitis resistance in a GWAS by using a high-density SNP array. Mastitis-resistant (15) and mastitis-susceptible (28) phenotypic extremes were identified from 224 lactating dairy cows on commercial dairy farm located in Utah based on multiple criteria of mastitis resistance over an 8-month period. Twenty-seven quantitative trait loci (QTLs) for mastitis resistance were identified based on 117 SNPs suggestive of genome-wide significance for mastitis resistance (p ≤ 1 × 10-4), including 10 novel QTLs. Seventeen QTLs overlapped previously reported QTLs of traits relevant to mastitis, including four QTLs for teat length. One QTL includes the RAS guanyl-releasing protein 1 gene (RASGRP1), a candidate gene for mastitis resistance. This GWAS identifies 117 candidate SNPs and 27 QTLs for mastitis resistance using a selective genotyping approach, including 10 novel QTLs. Based on overlap with previously identified QTLs, teat length appears to be an important trait in mastitis resistance. RASGRP1, overlapped by one QTL, is a candidate gene for mastitis resistance.

Metagenomic sequencing provides insights into microbial detoxification in the guts of small mammalian herbivores (Neotoma spp.).

Microbial detoxification of plant toxins influences the use of plants as food sources by herbivores. Stephen's woodrats (Neotoma stephensi) specialize on juniper, which is defended by oxalate, phenolics and monoterpenes, while closely related N. albigula specialize on cactus, which only contains oxalate. Woodrats maintain two gut chambers harboring dense microbial communities: a foregut chamber proximal to the major site of toxin absorption, and a cecal chamber in their hindgut. We performed several experiments to investigate the location and nature of microbial detoxification in the woodrat gut. First, we measured toxin concentrations across gut chambers of N. stephensi. Compared to food material, oxalate concentrations were immediately lower in the foregut, while concentrations of terpenes remained high in the foregut, and were lowest in the cecal chamber. We conducted metagenomic sequencing of the foregut chambers of both woodrat species and cecal chambers of N. stephensi to compare microbial functions. We found that most genes associated with detoxification were more abundant in the cecal chambers of N. stephensi. However, some genes associated with degradation of oxalate and phenolic compounds were more abundant in the foregut chambers. Thus, microbial detoxification may take place in various chambers depending on the class of chemical compound.

Long-range genomic regulators of THBS1 and LTBP4 modify disease severity in duchenne muscular dystrophy.

OBJECTIVE: Duchenne muscular dystrophy (DMD) is a severe X-linked recessive disease caused by loss-of-function dystrophin (DMD) mutations in boys, who typically suffer loss of ambulation by age 12. Previously, we reported that coding variants in latent transforming growth factor beta (TGFß)-binding protein 4 (LTBP4) were associated with reduced TGFß signaling and prolonged ambulation (p = 1.0 × 10-3 ) in DMD patients; this result was subsequently replicated by other groups. In this study, we evaluated whether additional DMD modifier genes are observed using whole-genome association in the original cohort. METHODS: We performed a genome-wide association study (GWAS) for single-nucleotide polymorphisms (SNPs) influencing loss of ambulation (LOA) in the same cohort of 253 DMD patients used to detect the candidate association with LTBP4 coding variants. Gene expression and chromatin interaction databases were used to fine-map association signals above the threshold for genome-wide significance. RESULTS: Despite the small sample size, two loci associated with prolonged ambulation met genome-wide significance and were tagged by rs2725797 (chr15, p = 6.6 × 10-9 ) and rs710160 (chr19, p = 4.7 × 10-8 ). Gene expression and chromatin interaction data indicated that the latter SNP tags regulatory variants of LTBP4, whereas the former SNP tags regulatory variants of thrombospondin-1 (THBS1): an activator of TGFß signaling by direct binding to LTBP4 and an inhibitor of proangiogenic nitric oxide signaling. INTERPRETATION: Together with previous evidence implicating LTBP4, the THBS1 modifier locus emphasizes the role that common regulatory variants in gene interaction networks can play in mitigating disease progression in muscular dystrophy. Ann Neurol 2018;84:234-245.